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It is shown that the p-wave singlet superconductivity, whose gap function is odd both 
in momentum and in Matsubara frequency, is reahzed prevaihng over the d-wave singlet 
state very close to the antiferromagnetic quantum critical point (QCP) both in the 
paramagnetic and in antiferromagnetic sides. Off the QCP in the paramagnetic phase, 
however, the d-wave singlet superconductivity with line-nodes is realized as conventional 
anisotropic superconductivity. This p-wave singlet superconductivity is an essentially 
gapless state, i.e., there is no gap in the quasiparticle spectrum everywhere on the 
Fermi surface due to its odd-frequency gap. These features can give a qualitative but 
nice understanding of the anomalous behaviors of NQR relaxation rate on CeCu2Si2 or 
CeRhlns where the antiferromagnetism and superconductivity coexist on a microscopic 
level. 

KEYWORDS: p-wave singlet superconductivity, gapless superconductivity, coexistence of antiferromag- 
netism and superconductivity, CeCu2Si2, CeRhlns 



1. Introduction 

The discovery of the superconductivity in CeCu2Si2^^ two decades ago triggered off a break- 
through in the field of superconductivity. Nowadays, we seem to have reached the consensus that 
the anisotropic spin singlet pairing, e.g. "d-wave singlet" pairing, is realized in the superconducting 
(SC) states close to the antiferromagnetic (AF) states due to the AF spin fluctuations.^"^) Recent 
developments on CeCu2Si2, however, suggest that it still disguises much various physics.^^^'''^^'^^'^^^ 
Polycrystalline sample of Ceo.99Cu2.o2Si2 shows a peak in the nuclear spin lattice relaxation rate 
1/Ti at SC transition temperature ~ 0.65K and 1/Ti oc T behavior at T < Tc, which indicates 
essentially gapless superconductivity to be realized at T < Tc at ambient pressure. At pressures 
P ^ O.lGPa, however, the 1/Ti oc behavior, indicating line-node gap, is observed. ^"^^ Since 
the Ge substitution for Si is considered to expand the lattice constant, giving negative pressure, 
1% Ge substituted compound CeCu2(Sio.99Geo.oi)2 exhibits AF ordering at T < Tn = 0.75K, fol- 
lowed by the onset of superconductivity at T = Tc = 0.5K at ambient pressure. The observation 
of single component of NQR signal, showing an appearance of the internal field throughout the 
sample, excludes the possibility of phase segregation between the SC and the AF phases; namely, 
SC and AF coexist on a microscopic level. In this SC states, 1/Ti does not show any significant 
reduction below but shows essentially the same behavior as the normal Fermi liquid state. At 
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Fig. 1. Schematic phase diagram near the phase boundary between antiferromagnetic and superconducting states 
based on the NQR experiments on CeCu2Si2 and CeRhlns. 

P = 0.85GPa, on the other hand, it does not exhibit the AF ordering and the hne-node SC gap 
evidenced from a clear behavior of 1/Ti oc ?^ 

Interestingly, the recently discovered pressure induced superconductor, CeRhlns of Ce-based 
heavy fermion, shows features quite similar to those of CeCu2Si2, while its crystal structure is 
rather two-dimensional compared to CeCu2Si2. At ambient pressure, CeRhlns shows AF ordering 
at T < Tn = 3.8 K and Tn exhibits a moderate variation for P ^ 1.75GPa, while the internal 
field Hi^t is linearly reduced in the region < P < 1.23GPa, suggesting that Hi^t vanishes at 
P ~ 1.6GPa. At P ~ 1.6GPa, the SC state begins to appear, even though the AF ordering still 
exists. ^^"^^^ The detailed NQR measurement shows that at P = 1.75GPa, the AF order and the SC 
order coexist on a microscopic level and 1 /Ti in SC state exhibits essentially the same T-dependence 
as the normal state as is seen in CeCu2Si2.^^-^^) Under the pressures P > 2.1GPa, where the AF 
state disappears, the 1/Ti oc behavior is observed as a usual singlet superconductivity with line- 
nodes. ^^-^ The specific heat of CeRhlns^^^ exhibits a feature consistent with the NQR measurement. 
At P ~ 1.6GPa, the specific heat shows a broad peak which is expected to be associated with the 
AF ordering at T ~ 3K, and shows a shoulder at T ~ 2K, which indicates the appearance of the 
gapless superconductivity. This shoulder grows to a mean- field like peak at P > 1.9GPa, which 
would be associated with the anisotropic singlet pairing. 

Putting all these together, we can draw a schematic phase diagram near the phase boundary 
between AF and SC states as shown in Fig.l. The P-T phase diagram can be separated into 
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three parts: AF phase, SC phase and a coexisting phase of AF and SC ordering (AF+SC phase). 
This coexistence is quite different from that of UPd2Al3, in which plural / electrons, in (5/)^- 
configuration of ion, shows a dual character of correlated electrons. Of three /-electrons per 
U"^^ ion, two electrons form essentially localized 5/^-electron state, as in the U^"*" state with singlet 
crystalline electric field ground state, while the remaining one electron forms the heavy Fermi liquid 
due to a larger hybridization with conduction electrons.^°^ In Ce-based heavy fermion compounds, 
on the other hand, there is only one 4/-electron per Ce^+ ion, and thus the same /-electron exhibits 
simultaneously itinerant and localized dual nature. The qualitative difference of the manner of 
the coexistence of AF and SC can easily be seen by the property of l/Ti on UPd2Al3^^) and 
CeCuaSia.^'^''^) 

A crucial aspect of the AF-I-SC phase is that the gapless superconductivity is realized, although 
the singlet pairing with line-node due to the antiferromagnetic fluctuation is expected to emerge. 
Although one can soon suspect that the origin of this gapless superconductivity is due to the 
impurity scattering, all the samples above are expected to be very clean. In fact, the clear l/Ti oc 
behavior is observed on the same sample under the pressures away from the boundary between 
AF-I-SC and SC phases. Even if the sample is very clean, there is still a possibility of the enhanced 
impurity scattering associated with quantum critical phenomena. ^^"^^^ In such a case, however, Tc 
should be reduced considerably, while Ceo.99Cu2.o2Si2 at ambient pressure, exhibiting 1/Ti oc T, 
does not show a clear decrease of T^.^^^ Therefore, the AF-I-SC phase in these Ce-based compounds 
seem to show the possibility of novel mechanism of "unconventional" superconductivity. 

In this paper, we present a possible scenario to unravel these mysterious features observed in 
CeCu2Si2 and CeRhIn5 near the boundary between AF-I-SC and SC phases. An essential idea here 
is that a gapless p-wave singlet superconductivity with the so-called odd-frequency gap, which is 
odd both in momentum and in frequency, should be realized very near the quantum critical point 
(QCP) and/or in the AF-I-SC phase. It depends on the specific property of the systems wheather 
the boundary between AF-I-SC phase and SC phase corresponds to QCP or not. However, the 
present theory can apply to both cases. The Pauli principle requires that the spin singlet gap 
function must be even under simultaneous parity- and time-inversion operations. Therefore, there 
are two types of singlet gap function; the first one is even both in frequency and in momentum, the 
other is odd both in frequency and in momentum. The latter type of gap was proposed by Balatsky 
and Abrahams (BA) in the context of the research of high-Tc cuprates about a decade ago,^^) after 
the proposal by Berezinskii for the model of superfluid ^He.^^^ BA argued that there is no gap in 
the quasiparticle spectrum with odd-frequency gap and these superconductors would still exhibit 
the Meissner effect, just as in the case of odd-energy gap proposed at very early stage of research as 
a model of CeCu2Si2 which was claimed to exhibit gapless nature of l/Ti.^^'^^) The model phonon 
interaction discussed by BA was shown to be an impractical one, and it was speculated without 
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specifying any explicit models that electron-electron interaction might mediate such odd frequency 
pairing in general.^^^ 

Here we show that a realistic model interction mediated by the critical spin fluctuatioin can give 
rise to the odd-frequency gap superconductivity if the Fermi surface of quasi particles satisfies some 
mild condition. In the next section, we describe the model and formalism of solving a linearized 
gap equation to obtain the SC transition temperature T^. In §3, we present the numerical results 
and discuss the condition that Tc of p-wave singlet superconductivity prevails over that of d-wave 
singlet superconductivity. In §4, we show that the obtained p-wave singlet SC state is really a 
gapless state also in real frequency space, i.e., A(a; = 0) = 0. In §5, we apply the same method to 
the region where the AF and SC order coexist and show the p-wave singlet superconductivity is also 
realized there. Then, we discuss the correspondence between results of our theory and experimental 
results. Section 6 is devoted to sammary. 



2. Model and Formulation 

First wc introduce the effective interaction taking into account the characteristic features of 
Ce-based compounds CeCu2Si2 and CeRhlns near the boundary between the AF and SC phase. 
For this purpose, it may be useful to recall the fact that the SC transition temperature of CeCu2Si2 
(CeRhlns) exhibits maximum at away from the phase boundary P* ~ 2.5GPa^^'^^) (2.5GPa)). 

It was shown theoretically, on the standard approximation for an extended Anderson lattice 
model, that the maximum of Tc is raised by the critical valence fluctuations. The recent exper- 
iment on CeCu2(Sio.9Geo.i)2 has revealed that the SC phase is separated into the phase near the 
phase boundary of AF phase and that around P ~ P* }'^) So, it is natural to consider that the latter 
is due to the critical valence fluctuations and the former is due to the critical spin fluctuations. This 
would be also valid for the CeRhlns which exhibits the analogous features with the CeCu2Si2.^^^ 
Upon this consideration, in order to discuss the mechanism of superconductivity near the phase 
boundary, we introduce the effective interaction ^(q, iwm) the following phenomenological form 
mediated by spin fluctuations: 

ri + Aq +C\LOm\ 

where g is the coupling constant, Np the densiyy of states at the Fermi level, and = 4 -|- 
2(cosq'j; -I- cosqy) in two dimensions. This type of pairing interaction was adopted by Monthoux 
and Lonzarich to discuss the strong coupling eflFect on the superconducticity induced by the critical 
AF fluctuations.^) One of the motivations of adopting q^ as above is that the ordering vector 
Q of AF phase in CeRhlns was identified as Q = (1/2,1/2,0.297) by the neutron diffraction 
measurements,^^) while that in CeCu2(Sio.99Geo.oi)2 has not been identified yet. Another one is 
that calculations become much simpler in two dimensions than three dimensions while qualitative 
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physical picture is remained unchanged. 

The linearized gap equation is given in a weak-coupling approximation as follows: 

A(k,is„) = -T ^ ^^^:2'':f:i;''"^ A(kM4), (2) 

where is the dispersion of the quasiparticle, and tlie pairing interaction is given by (1). The 
pairing interaction can be decomposed as 

V{k - k', iu;m) = ^ Vi(ic^„)</.r(k)<^Kk'), (3) 

I 

where (pi(\i) is an irreducible representation of crystal group. The coefficient Vi{iuJra) in (3) is given 
as 

ViiiUm) = J]</.Kk)y(k - k',iuJm)4>li^')- (4) 
k,k' 

In order to make our argument as simple and transparent as possible, we concentrate on the 
frequency dependence of the gap function and assuming the two-dimensional circular Fermi surface 
shown in Fig. 2. Then, the p- and d-wave components are retained: 

/r/k Jk' 
</)^(k)F(k-k',ia;^)<^r(k), (5) 
VF VF 

(^p(k) = ^-^N-H{^i, - n) sin kFx {= z-^N-H{i], - n) sin kFy), (6) 
0d(k) = $;^^iV^^(5(Ck - M)(cosfeFx - cos/cFy), (7) 

where the wave vector on the Fermi surface are expressed as k^? = {kFx ^ kFy), and is the 
normalization factor. 

The linearized gap equation for each partial-wave component can be written in the form: 

A(r)AKis.) = -tY: (8) 

Assuming = "^^dkl — kp), one can carry out the integration of eq. (8) with respect to k' as 
follows: 

1 rkp+ko 

dkp- 



'^'^ Jkp-ko v%{k-kFf + \en\^ 

kF . _i VFko 

tan - — -. (9) 



'KVF\en\ \£n\ 

As mentioned above, the gap function must satisfy the following symmetry for the spin-singlet 
pairing 

A<i(k, lEn) = Arf(-k, iSn) = Ad(k, -iSn) = Ad(-k, -ie„), (10) 

Ap(k,i£„) = -Ap(-k,i£„) = -Ap(k,-i£„) = Ap(-k,-i£„). (11) 
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Fig. 2. Relation between the Fermi surface (kp = 0.75-k) and the ordering vector Q. 

Here we measure energies in unit of the Fermi energy ep. g^Np is set to be equal to ep throughout 
the paper. The Fermi surface (FS) is assumed to be two-dimensional circular with kp = O.TSvr 
as is shown in Fig. 2. This simplified model may give clear understanding of a condition for the 
realization of the odd-frequency gap, although it does not corresponds to the practical Fermi surface 
of the heavy fermion compounds. 

The transition temperature Tc is determined by condition A(Tc) = 1. We have solved the eigen- 
value problem (8) numerically by retainning 512 Matsubara frequency up to |en| < 1315837rT. 
The temperature dependences of the eigen value A(T)'s are shown in Fig. 3. For the d-wave pair- 
ing, A(T) increases monotonously (logarithmically) as the temperature decreases as usual. For the 
p-wave odd-frequency pairing, on the other hand, A(T) exhibits a broad peak as shown in Fig. 3. 
Since this feature can be seen even in the case of the much larger cut-off frequency, it is not due to 
the cut-off frequency, but due to an intrinsic nature of the odd-pairing in frequency. Consequently, 
it is possible that there exist two temperatures satisfing X{T^) = 1, so that the p-wsve pairing 
is realized only in the temperature region T^T < T < (see the case r] = 0.05 in Fig. 3). But in 
the case tj < 0.04, A does not fall below unity even for T — > 0, so that the p-wave pairing is realized 
in < T < r+. 

We recall that BA reached the anologous results for the odd-frequency pairing that the normal 
phase reappears below the lower-Tc. In the present model, we cannot show the existence of T^, 
since the transition temperature for d-wave singlet pairing is higher than for r] > 0.04. 
The detailed analysis about the possibility of this multiple transition will be given in a subsequent 
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Fig. 3. Temperature dependence of the eigen value A for p-wave singlet (p-S) and d-wave singlet (d-S). 
paper. 

3. Competition between p-wave and d-wave pair states 

As is shown in Fig. 4, the transition temperature for p-wave singlet Tc exceeds when 
r] ^ 0.02. That is to say, in the region away from the QCP, d-wave singlet pairing is realized and 
p-wave singlet pairing is formed very close to the QCP. 

The reason why such p-wave singlet pairing arises near the AF-QCP can be understood as 
follows. The frequency dependence of Vi{iuJm), which is calculated numerically by eq. (5), are 
shown in Fig. 5. As is also shown in Fig. 5, it should be reasonable to approximate the frequency 
dependence of V/(iWm) as 

Vi{iuJm) ^ VI In -^-^ — -. (12) 

Considering the symmetry property of the gap function, eqs. (10) and (11), we can express each 
Vi in the form: 




(13) 
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Let US consider the most singular case, £„ = = ttT. Of course this example does not 
reproduce the present situation exactly, but gives a simple picture of the competition of the pairing 
states. Each Vi is given by 

Vi(0) = fln^, (15) 
When the system locates away from the QCP, i.e., fj S> 2TrT, 



Vp2TrT 




(17) 

_21n?7-^). (18) 
V J 

In this case, Vp is almost negligeble compared to V^, leading to S> Tc- 

In contrast, when the system locates very close to the QCP, i.e., fj ^ 2TrT, 

- (19) 




Va{0) ~ _21n(27rr)-^-ln^U-^ln^. (20) 

'^^ ' 2 I ^ ' 2ttT 27rT J 2 2'kT ^ ' 

This time Vp is comparable to V^, and in the case of \vp\ > \vd\, 7? becomes higher than T^. 

As shown in Fig. 6, \vi\ considerably depends on kp- In other words, it is determined by 
the relation between the Fermi surface and AF ordering vector Q. In the present case, the most 
singular pair scattering with q = Q = (=b7r, ±7r) uses the points near (±kF/V2,±kF/V2), where 
the nodes of the d-wave pairing exists while that of the p-wave pairing does not. (See Fig. 7 
(a), (b)) Therefore the effects of the singular pair scattering are suppressed in the case of d-wave 
pairing. This causes the relation \vp\ > \vd,\ in the present case. 

Considering the above discussion, we can extend the condition of the emergence of the p-wave 
singlet pairing to the general case as following. Firstly, the pair scattering interaction must have a 
sharp peak in frequency whose width is smaller than the temperature (in the present case, r] <^T). 
Secondly, the most singular pair scattering with vector Q has to be canceled out by the nodes 
of the d-wave pairing. In the case of Hubbard model at the half-filling, for example, the Fermi 
surface is nested with Q = (±7r, itTr) and the suppression of pair scattering becomes small, namely, 
\vp\ ~ l^^dl (see Fig 7 (c)). So, it would be difficult for the p-wave singlet pairing to prevail over 
the d-wave singlet pairing on the Hubbard model, even near the AF-QCP. How about the cases of 
CeCu2Si2 and CeRhlns? The l/Ti of both compounds rapidly increase approaching the AF+SC 
phase from SC phase. ^' This indicates that the system is located very near the QCP, which would 
satisfy the first condition mentioned above. In CeRhIn5, the most singular pair scattering vector is 
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Fig. 6. A;F-dependence of vi. The vertical dotted-line indicates kp = 7r/\/2 where |Q| is equivalent to the diameter 
of the Fermi surface. 
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Fig. 7. Illustration of the relation between the AF ordering vector Q and the position of the nodes, (a) p-wave 
pairing case: the most singlular pair scattering with textbfQ is not suppressed by the existence of the nodes, (b) 
d-wave pairing case: the scattering with Q is suppressed by the existence of the nodes, (c) d-wave pairing with FS 
nested with Q case: the suppression of pair scattering becomes small. 
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Fig. 8. Frequency dependence of the gap function of the p-wave singlet pairing at T+ /ep = 7.0 x 10"^ and /ef = 
1.0 X 10~'*for r] = 0.04. The sohd hne indicates the results of Fade approximation. 

Q = (vr, vr, 0.2977r).^^) Its Fermi surface does not nest in the direction Q and is satisfied the relation 
ki? — = Q in the direction (±1, ±1, 0).^^^ The feature that Q is independent from the nesting 
vector can be understood from the dual nature of electron. According to the itinerant-localized 
duality model, the spin susceptibility xiQ-, i'^m) is expressed in terms of the polarization function 
n(g, \ujm) of quasiparticles and the exchange interaction J(q, iLOm) between localized component of 
electrons as follows: 

X(g,iu;m)~"^ = Xo(i'^m)"^ - 2\^Il{q,\uim) - J{q,iuJm), (21) 

where xo is the local susceptibility expressing an effect of Kondo-like correlation, and A is the 
renormalized spin-fermion coupling. Even if the Fermi surface is not nested, i.e., 11 is not enhanced, 
the AF ordering can be triggerd by J(Q). Therefore, CeRhlns would satisfy both conditions. In 
the case of CeCu2Si2, a situation is much more complicated (Q ~ (0.277r, 0.277r, 0.527r)3^) ), but its 
Fermi surface seems to satisfy the condition. 

4. Gaplessness of odd-frequency pairing 

The frequency dependence of the gap function of the p-wave singlet pairing at = 7.0 x 10^^ 
and T~ = 1.0 x 10~^ for rj = 0.04 are shown in Fig. 8. For T = T^, at first sight, the gap seems 
to be divergent when e„ approaches zero. The gap at T = T~ , on the other hand, has a peak at 
e* = 6.19 X 10~^, which is comparable to T^, and vanishes for e ^ 0, i.e., it exhibits a gapless-like 
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structure. This feature seems robust against variations of 77, kp or g. Therefore, we can suspect 
that the peak and the gapless-hke structure is only hidden for T = T+. 

At the present stage, however, we cannot declare wheather the present p-wave singlet gap 
function is completely gapless or not. So, we have to perform the analytic continuation of the gap 
function. In order to do that, first we apply the Pade approximation to the calculated gap function. 
Since the gap function is expected to be proportional to |en|~^ for large £„, we take the form: 

1 + a2£„ + 036^ + 046^ 

for En > 0. We apply the least squares method to the calculated gap function A(i£„)by Ap(i£„) and 
obtain ai = -1.2455 x lO^^, 02 = 2.2966 x 10, 03 = 1.6980 x 10^, 04 = 7.7980 x 10^, 05 = 6.6046 x 10. 
Using these values of 01^5, we display Ap(i£„) in Fig. 8, in which Ap(i£„) shows a good agreement 
with A(i£„). The analytic continuation of this form becomes 

ReAp(£) = — ^ — , (23) 

(1 - a3£^)^ + £^(02 - a4£^) 

a5£(a3£2 - 1) + aia5e(a2 - a4£^) 
'"^^^^^ = (l-a3£^)^ + £^(a2-a4£^) " ^^'^ 

From eqs. (23) and (24), the gap function at £ = is reduced to ReAp(O) = —8.2260 x 

10~^, ImAp(O) = 0. Therefore, we can conclude the present p-wave singlet superconductivity is a 

technically gapless one. We can determine the maximum of the gap, Ao(T) below T^■. by assuming 

that the gap has the same frequency dependence as that of the eigenfunction at T ^ Tc, where 

the eigenvalue A exceeds unity further. Namely, A(i£„;T) = Ao(r)Ap(i£„; T <C Tc). At T=0, we 

obtained Ao(0)/Tc = 0.42 (by setting parameters as ai = -1.5256 X 10~^,a2 = 4.6835 x 10, as = 

6.4780 X 10^,0,4 = 2.8209 x 10^, = 1.1492 x 10^, which corresponds to r? = 0.01 case), which is 

much smaller than the BCS relation Ao(0)/Tc = 1.765. 

We can also discuss the gapless nature in more general form as follows. The linearized gap 

equation for the odd-frequency gap can be written as 

A,(i£„) = -Np Tdgr^ l[A,(i£;)- A,(-i£;,)] (25) 

-'00 , Pnl "T S ^ 

r°° dr T 1 

= -NpP / — coth [Vi-x - iS) - V{-x + iS)] X [A(i£„ + x) - A(-i£„ - x)] 

dx X 1 

The analytic continuation, i£„ — £ + id in eq. (26), results in the following gap equation: 

r'~^ dx ( X — E 

A(£ + i5) = -NfP j — |coth^^[F(£-x-i5)-F(£-x + i(^)] 

+ tanh ^ [V{e - x + \5) -V{e + x + \5)\ } 
A(x + i5) - A(-x - i5) 
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V{u + i6) = vi{lrL 



—1 \ 

, + 1 tan — 

Va;2 + ^2 r]^ 



(28) 



Substituting eq. (28) into eq. (27), we obtain 



(e + x)2 + 



(e — x)2 + 



+1 2 coth — — — tan h tanh — - ( tan tan 

A(,r + if)) - A(-.r - i5) 



V 



)]} 



(29) 

Let us express the real (imaginary) part of the gap as A'(£) (A"(£)), i.e., A(£ + id) = A'{e) + 
iA"(e). Then, each part of (29) are reduced to 



A'(s) = - 



2 coth 



e — X 
2T 



tan 



e — X 



(£ + x)2 + f)2 A'{x)-A\-x) 
{e — x)2 + 1^2 x 

X , _i £ — X 



+ tanh — (tan ^ 



»7 



. _ie + x. 
tan — r — j 



A"(x) + A"(-x) 



(30) 



(£ + x)2 + r?2 A^^(x) + l^"{-x) 
{e — x)2 + 7^2 



£ — X 

2 coth — — — tan" 



2T ri 
A'(x) - A'(-x) 



A S X ^ X / 1 

+ tanh — (^tan 



X 

S X 1 ^ ~l~ ^ \ 

— z tan — z — j 



X 



+ {A"(x) + A"(-x)}7r5(a 



(31) 



Note here that we can show easily A'(£) = -A'(-£), A(£) = A"(-£). At £ = 0, 

A'(0) = 0, (32) 
A"(0) = -viN^P P ^ (^coth ^ - tanh ^) tan^i | • ^ - viNp^^/^" {Q). (33) 

Prom eq. (33), 

T\ ^n,^^ f^dx/ , X , x\ ix A'(x) 

1 + viNftt^] A"(0) = -viNpP J — (coth — - tanh — J tan'^ -z ■ (34) 



The right hand side of eq. (34) is constant, and 1 + viNp-KTjr] — — oo when 77 — 0. Therefore, in 
order to satisfy eq. (34), 



A"(0) ^0 for 77 ^ 



(35) 



is needed. Prom eqs. (32) and (35), we reach the fact that the present p-wave singlet gap function 
leads to an essentially gapless superconductivity at the QCP. 
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Considering the low-frequency structure of the gap function, we assume that the gap function 
can be approximate in the form: 

A(k,iu;^) = ^ia;^(/!>p(k). (36) 
Then we find from the poles of Green's function that 

= , ^ ■ (37) 



^l + (Ao/re)202(k) 

Therefore, the quasiparticle spectra with such gap function is gapless, i.e., there is no difference 
in excitations between the nomal states and the present SC states, except for the effective mass 
enhancement: 

ml = m^l + (Ao/re)202(k). (38) 

This mass enhancement is, however, very small due to the smallness of Aq/Tc ~ 0.42 as mentioned 
above, so that the specific heat or the NMR relaxation rate does not show any significant change 
at Tc. This would correspond to the l/Ti behavior observed in Ceo.99Cu2.o2Si2^~^^ which shows 
almost the same behavior as the normal Fermi liquid state. If the 1/Ti oc T behavior is due to 
the impurity scattering, l/Ti should shows a significant reduction at and exhibit l/Ti oc T well 
below Tc.^^^ This is the reason why we conclude the l/Ti oc T is not due to the impurity scattering 
but due to the odd-frequency gap. 

5. Coexistence of AF and SC order 

In the previous sections, we discussed the emergence of the gapless p-wave singlet pairing 
prevailing the d-wave singlet pairing in the paramagnetic (PM) backgrounds. The experimental 
results^"^' 16-19) suggest that the gapless superconductivity is realized in the AF backgrounds rather 
than in the PM backgrounds. (In the SC phase very close to the phase boundary, the present theory 
predicts that the gapless ^>wave singlet superconductivity is realized, but the detailed experimental 
results in this region have not been obtained yet.) Here we discuss that the gapless p-wave singlet 
can also be realized in rather wide region in AF-I-SC phase. Since the transverse spin susceptibility 
X± would be much dominant than the longitudinal one xy , even if the damping effect would have 
been taken into account, the pairing interaction in the AF backgrounds may given as follows: 

2 g'^Np 

V{q, iujm) = g X±(q, it^m) = 2-2.1 fi' (39) 

where S corresponds to the spin-wave velocity. This interaction is similar to that just at QCP. 
Namely, the frequency dependence of ^(^^(ia;^) can be approximated as 

Vf^^{iu:m)^vi\n-^. (40) 
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Fig. 9. Illustration of the frequency dependence of Vi{iujm)in the AF backgrounds, where lj^ = e„ ± e^. The p-wave 
channel dominates the d-wave channel because Vp > Vd and the unsaturated behavior of Vi{iLj,ri)- 

Comparing this with the expression (12), the interaction V/^^ of the AF side can be obtain by 
the transformation Vi^^{iuJm) 2Vj^^ {iu! m','fi = 0). Figure 9 illustrates this situation. There is 
no saturation in contrast to V^^'^ in the PM side (see Fig. 5 ), which stimulates the emergence 
of p-wave singlet pairing. How the p-wave states dominate the d-wave one can be understood in 
the similar scenario discussed in §3, but this time we present much more generalized and intuitive 
picture about how the p-wave dominates the d-wave. Each Vj^^ can be expressed in the form 



where uj^ = En ^ £n- From these expressions, we can easily find that V^^ > Vj^^ , in general, for 
Vd ~ Vp because the first term and the second term are added for d-wave, while they offset each 
other for p-wave. In the case Vp is moderately larger than Vd as displayed in Fig. 9, however, 
the first term of V^^^ is much larger than V^^, so that the p-wave singlet pairing dominates the 
d-wave one. This situation would be realized easier in the AF side rather than that in the PM side 
because of the factor 2 and the unsaturation behavior of V/^^ . Judging from above, no matter how 
the detail is, the p-wave pairing prevails d-wave one when the frequency dependence of the pairing 
interaction takes the form as shown in Fig. 9. 

Below Tn, however, the energy gap is formed around the hot spot (Fig. 10) and the missing 
part of the Fermi surface spreads out as the AF ordering is developed. It leads to the remarkable 




(41) 



(42) 
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hot spot 




Fig. 10. Missing Fermi surface due to the AF energy gap and the definition of (p. 

reduction of because the dominant pair scattering is suppressed due to the existence of the AF 
gap. To simulate this situation, we parameterize the effect of the AF order as the angle (/? of the 
missing part of the Fermi surface as displayed in Fig. 10. The growth of the AF gap, spreading of 
deprives the pairing interaction of the most singular scattering with the AF ordering vector Q, 
and suppresses V/^^ . It is noted that this suppression is not the saturation effect displayed in Fig. 
5, but the reduction of vi. This reduction of vi is remarkable for fp, because the enhancement of 
Vp is mainly due to the scattering using the hot spots. 

We can also discuss the gapless nature as discussed in §4. In the AF backgrounds, the pairing 
interaction V^^^ is given by the same form as that in the PM backgrounds with vi — > 2vi in the 
limit r/ ^ 0. Thus we reach the same conclution that the gap function in the AF side is always 
gapless: 

AkF(O) = 0, (43) 
Aaf(O) = 0. (44) 

Together with the T^s in the PM phase, the calculated T^s near the QCP are shown in Fig. 
12 Here we set S = 1.2 and the other parameter g^Np and kp to be the same as those in the 
PM backgrounds. In the PM region, the distance from the QCP is parameterized by ??, and 
in the AF region, that is parameterized by ip. Indeed, the gap in AF backgrounds, shown in 
Fig. 11, has a gapless-like structure, and the parameters of Pade approximation are obtained as 
ai = -3.3415 X 10^3,02 = 2.5624 x 10,03 = -7.7475 x 10,04 = 3.3734 x 10^, 05 = 4.8612 x 10. 
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Fig. 11. Frequency dependence of the gap function of the p-wave singlet pairing in the AF background. The soUd 
hne indicates the resuhs of Fade approximation. 
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Fig. 12. Transition temperatures near the QCP. AF+dSC: d-wave singlet pairing (d-S) with hne nodes; AF+pSC: 
gapless p-wave singlet (p-S); pSC: gapless p-S; dSC: d-S with hne nodes. 



18 Yuki FuSEYA, Hiroshi KoHNO and Kazumasa Miyake 

The gapless odd-frequency p-wave singlet pairing is realized only near the QCP, and d-wave 
singlet pairing prevails p-wave singlet pairing in the PM region away from the QCP. According to 
this phase diagram, the observed l/Ti may be explained as follows. Near the QCP, i.e., AF+pSC 
and pSC region displayed in Fig. 12, l/Ti does not show any significant reduction due to the gap- 
less SC state. Away from the QCP (AF+dSC and dSC region), l/Ti oc behavior is observed. 
This phase diagram agrees qualitatively with several experiments in CeCu2Si2 and CeRhlns. Poly- 
crystalline sample of Ceo.99Cu2.o2Si2 shows gapless behavior at ambient pressure, and l/Ti oc 
behavior at pressures P ^ O.lGPa.^'^^ This property can be understood that Ceo.99Cu2.o2Si2 at 
P = OGPa is located inpSC phase and at P ^ O.lGPa, dSC phase. CeCu2(Sio.99Geo.oi)2 compounds 
shows the AF ordering at T < 0.75K, and below 0.5K, the gapless superconductivity, coexisting 
with the AF order, at ambient pressure. This compounds also exhibits the line-node gap under 
the pressure P = 0.85GPa.^^ Therefore, this compounds would be located in AF+pSC phase at 
ambient pressure and in dSC phase at P = 0.85GPa. The existence of pSC phase in this sample 
cannot be recognized from the present experimental data. Similarly, in CeRhlns the pressure region 
1.6 < P < 1.75GPa would correspond to AF+pSC and P > 2.1GPa would correspond to dSC.^^^^) 
For each case, the AF+dSC phase, where the d-wave singlet SC states and the AF states coexist, 
have not been observed yet. 

6. Conclusion 

In the present paper, wc have shown that the p-wave singlet superconductivity with the gap 
function which is odd in both momentum and frequency prevails the d-wavc singlet superconduc- 
tivity, at the AF-QCP and in the AF side of AF-QCP. This odd-frequency p-wave singlet pairing is 
stabilized in the AF states, and is able to coexist with the AF order. The characteristic properties 
of this p-wave singlet superconductivity is that i) there is no gap in the quasiparticle spectrum; ii) 
d-wave pairing arises apart from the AF-QCP even in the AF state; iii) In some parameter region, 
p-wave pairing is realized only in the window of temperatures T~ < T < . The first and the 
second property would explain the experiments. Namely, the 1/Ti oc T behavior is observed in the 
AF+SC region or in the boundary of CcCu2Si2 or CeRhlns, and 1/Ti oc T"^ behavior in the PM 
side away from the boundary. From the second one, we can predict that the d-wave pairing will be 
observed again away from the AF-QCP in the AF states. 

The condition of the emergence of the p-wavc singlet pairing is that the FS is not nested, so 
that the hot points of surface connected by the AF wave vector Q = (vr, vr) are isolated, while 
the AF order is induced by the exchange interaction between localized component of spin degrees 
of freedom. Although the structure of the Fermi surface of CcCu2Si2'^^^ and CcRhlns'^"^^ is much 
more complicated than the present model, they seem to satisfy the condition from inspection of 
the shape of the Fermi surface obtained by band structure calculations. We thus believe that the 
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present theory would distill the fundamental picture of both CeCu2Si2 and CeRhIn5, and obtain 
the qualitative understanding of them. In order to get futher understanding, extensive calculation 
considering the momentum-dipendence and the practical dispersion of the quasiparticle are now in 
progress. 
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